Counter Terms for Low Momentum Nucleon-Nucleon Interactions 
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There is much current interest in treating low energy nuclear physics using the renormalization 
group (RG) and effective field theory (EFT). Inspired by this RG-EFT approach, we study a low- 
momentum nucleon-nucleon (NN) interaction, Vi ow -k, obtained by integrating out the fast modes 
down to the scale A ~ 2fm -1 . Since NN experiments can only determine the effective interaction 
in this low momentum region, our chief purpose is to find such an interaction for complex nuclei 
whose typical momenta lie below this scale. In this paper we find that Viow-h can be highly 
satisfactorily accounted for by the counter terms corresponding to a short range effective interaction. 
The coefficients C n of the power series expansion ^ C n q n for the counter terms have been accurately 
determined, and results derived from several meson-exchange NN interaction models are compared. 
The counter terms are found to be important only for the S, P and D partial waves. Scaling behavior 
of the counter terms is studied. Finally we discuss the use of these methods for computing shell 
model matrix elements. 

PACS numbers: 21.60.Cs; 21.30.fte; 27.80.+j 



I. INTRODUCTION 

Since the pioneering work of Weinberg [|, there has 
been much progress and interest in treating low-energy 
nuclear physics using the renormalization grpur 
and effective field theory (EFT) approach 0, y 
H II H H 13 113 A central idea here is that physics 
in the infra red region must be insensitive to the details 
of the short range (high momentum) dynamics. In low- 
energy nuclear physics, we are probing nuclear systems 
with low-energy probes of wave length A; such probes 
cannot reveal the short range details at distances much 
smaller than A. Furthermore, our understanding about 
the short range dynamics is still preliminary, and model 
dependent. Because of these considerations, a central 
step in the RG-EFT approach is to divide the fields into 
two categories: slow fields and fast fields, separated by a 
chiral symmetry breaking scale A x ~ 1 GeV. Then by in- 
tegrating out the fast fields, one obtains an effective field 
theory for the slow fields only. This approach has been 
very successful in treating low-energy nuclear systems, as 
discussed in the references cited above. 

In order to have an effective interaction appropriate 
for complex nuclei in which typical nucleon momenta are 
< hp, we separate fast and slow modes by a scale much 
smaller than A x , namely A ~ 2 fm . This corresponds 
to the pion production threshold, a scale often not con- 
sidered in traditional EFT methods, even though it rep- 
resents the limit of available 2-nucleon data; experiments 
give a unique effective interaction only up to this scale. 
This cutoff has been e mployed recently by several authors 
dill, 13 13 13 III 13 III HI in a new development 
for the nucleon-nucleon (NN) interaction: the derivation 



of a low momentum NN potential, Vi ow -k- Despite some 
general similarities, we would like to stress that Vi ow -k is 
not a traditional EFT construct, and would like to now 
attempt to clarify this issue. 

One drawback to the commonly used RG-EFT, which 
adheres to order-by-order consistency in the power count- 
ing, is it's lack of predictability. The number of unknown 
parameters increases rapidly as the number of nucleons 
involved in the process increases. Even if such a method 
were available for complex nuclei, the parameters would 
be far from first principle QCD. The key ingredient in our 
approach is the highly refined standard nuclear physics 
approach (SNPA) and the objective is to marry the SNPA 
to an EFT |3|. 

As an illustration of our strategy, we briefly discuss 
how this marriage can be effectuated. As summarized 
recently [3 |2^, a thesis now develped for some time 
posits that by combining the SNPA, based on poten- 
tials fit to experiments, with modern effective field the- 
ory, one can achieve more predictive power than the 
purist's EFT alone. This combination of SNPA and ef- 
fective field theory, called EFT* by Kubodera (H, but 
more properly called more effective effective field theory 
(MEEFT) as noted by Kubodera, was rediscovered in 
[3 13 El 13 13 13 111 I3H3, and consists of anRG- 
type approach which limits A to 2.1fm~ 1 , because this is 
equal to the center of mass momentum up to which exper- 
iments measuring the nucleon-nucleon scattering phase 
shifts have been carried out. 

Thus, the "more effective" in MEEFT arises from the 
guarantee that all experimental data will be reproduced 
by Vi ow -k, which will contain no Fourier components 
higher than those at which experiments have been car- 
ried out and analyzed. Therefore, Vi ow -k is reliable up 
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to k = A. This seems to be adequate for a description of 
shell model properties of complex nuclei. 

In practice, the construction of Vi ow -k begins with 
one of a number of available realistic models 0, |H 
0111 El for the nucleon-nucleon (NN) potential V NN . 
While these models agree well in the low momentum 
(long range) region, where they are just given by the one 
pion exchange interaction, in the high momentum (short 
range) region they are rather uncertain and, in fact, dif- 
fer significantly from one another. Naturally, it would 
be desirable to remove these uncertainties and model de- 
pendence from the high momentum components of the 
various modern NN potentials. Following our RG-type 
approach to achieve this end, it would seem to be ap- 
propriate to integrate out the high momentum modes of 
the various Vnn models. This is how the Vi ow -k was de- 
rived, and in section II, we provide a brief outline of the 
derivation based on T-matrix equivalence. 

The low-momentum NN potential, Vi ow -k, reproduces 
the deuteron binding energy, low-energy NN phase shifts 
and the low-momentum half-on-shell T-matrix. Further- 
more, it is a smooth potential and can be used directly 
in nuclear many body calculations, avoiding the calcula- 
tion of the Brueckner G-matrix El El El El El HI 
Shell model nuclear structure calculations u sing Vipw-h 
have indeed yielded very encouraging results |13L fl9L l20j . 
These calculations apply the same Vi ow -k to a wide range 
of nuclei, including those in the sd-shell, tin region, and 
heavy nuclei in the lead region. At the present, no alge- 
braic form exists for Viow-k) an d it should be useful to 
have one. In the present work we study the feasibility of 
finding such an expression. 

A central result of modern renormalization theory is 
that a general RG decimation generates an infinite series 
of counter terms H 0, H HH consistent with the input 
interaction. When we derive our low momentum interac- 
tion, the high momentum modes of the input interaction 
are integrated out. Does this decimation also generate a 
series of counter terms? If so, then what are the proper- 
ties of the counter terms so generated? We study these 
questions in section III where we carry out an accurate 
determination of the counter terms and show that this 
approach reproduces not only Vi ow -k, but the deuteron 
binding energy and low-energy phase shifts as well. We 
shall discuss that the counter terms represent generally a 
short range effective interaction and are important only 
for partial waves with angular momentum I < 2. The 
scaling behavior of the counter terms with respect to the 
decimation momentum will be studied in section IV. Fi- 
nally in section V, we examine the prospect for using 
Viow-k and the counter term method in shell model cal- 
culations. 



II. T-MATRIX EQUIVALENCE 

Since our method for deriving the low-momentum in- 
teraction Viow-k has been described elsewhere 0, 0, 
IbH ]. in the following we only outline the derivation. We 
obtain Viow-k from a realistic Vnn model, such as the 
CD-Bonn model, by integrating out the high-momentum 
components. This integration is carried out with the 
requirement that the deuteron binding energy and low- 
energy phase shifts of Vnn are preserved by Vi ow -k- This 
preservation may be satisfied by t he following T-matrix 
equivalence approach 0, 0, [l4| . 
half-on-shell T-matrix 



We start from the 



T(k', k, k 2 ) = V NN (k', k) + / q 2 dqV N N(k', q) 

Jo 

X T2 T- ^T(q,k,k 2 ), (1) 

k z — q z + z0 + 

noting that the intermediate state momentum q is inte- 
grated from to oo. We then define an effective low- 
momentum T-matrix by 



Tiow-k{p',P,P ) = Vi ow - k (p',p) + / q dqVi ow -k(p',q) 

Jo 

x ~ 2 , . n+ Tiow-k(q,P,p 2 ) (2) 

p — q + «U + 

where A denotes a momentum space cut-off and (p',p) < 
A. We choose A ~ 2fm _1 , essentially the momentum 
up to which the experiments give us information in the 
phase shift analysis. Note that in Eq.(2) the intermediate 
state momentum is integrated from to A. We require 
the above T-matriccs satisfying the condition 

T(p',p,p 2 ) = Ti ow ^ k (p',p,p 2 ); (p',p) < A. (3) 

The above equations define the effective low momentum 
interaction Vi ow -k, and are satisfied by the solution 0, 

mm 



v l0 



Q' / Q 



Q / Q 



-Q' Q Q Q 



(4) 



which is just the Kuo-Lee-Ratcliff folded-diagram effec- 
tive interaction j3l l3l| . Here Q represents the irre- 
ducible vertex function whose intermediate states are all 
beyond A; Q' is the same as Q except with its terms first 
order in the interaction removed. 

For any decimation momentum A, the above Viom-k 
can be calculated highly accurately (essentially exactly) 
using cither the Andreozzi-Lee-Suzuki 

mm (als) or 

the Krenciglowa-Kuo |3J| iteration methods. These pro- 
cedures preserve the deuteron binding energy in addition 
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FIG. 1: Comparison of Vi ow -k with PVnnP plus counter 
terms, for 1 S'o and 3 Si channels. 

to the half-on-shell T matrix, which implies preservation 
of the phase shifts. After obtaining the above Vi ow -k> 
which is not Hermitian, we further perform an Okubo 
transformation .353 to make it Hermitian. This is done 
using the method given in Refs. j3|||37|]. Here we first 
calculate the eigenvalues and eigenvectors of the opera- 
tor oj^u) where u> is the wave operator obtained with the 
ALS method. Then the Hermitian Vi ow ~k is calculated in 
terms of these quantities in a convenient way (see Eq.(23) 
of Ref.|13). 

III. COUNTER TERMS 

Here we study whether the low-momentum interaction 
Viow~k can be well represented by the low-momentum 
part of the original NN interaction supplemented by cer- 
tain simple counter terms. Specifically, we consider 

Viow-k(q, q') — Vbare{q, q') + V coun ter(q, q'); (<?, q') < A 

where Vbare is some bare NN potential Vnn- Since 
Viow-k and Vbare are both known, there is no question 
about the existence of V coun ter in general. Our aim is, 
however, to investigate if V counter can be well represented 
by a short ranged effective interaction, such as a smeared 
out delta function. We shall check this conjecture by as- 
suming a suitable momentum expansion form for V coun ter 
and investigating how well it can satisfy the above equal- 
ity. Since Vbare is generally given according to partial 
waves as is Vi ow -k, we shall proceed to determine V coun ter 
separately for each partial wave. At very small radial dis- 
tance, the Bessel function ji(qr) behaves like (qr) 1 . We 

aSSUme that Vcoimier 

is a very short ranged interaction. 




k(fm ) 

FIG. 2: Comparison of Vi ow -k with PVnnP plus counter 
terms, for 3 Si— 3 Di and 3 Di channels. 

Hence the leading term in a momentum expansion of the 
partial wave matrix element {ql\V coun ter\q'l'} is propor- 
tional to q l q a . Thus we consider the following expansion 
for the partial wave counter term potential, namely 

(ql\V counter \q'l') = q l q ,l '[C + C 2 (q 2 + q' 2 ) 

+C,{q 2 + q' 2 ) 2 + C^ + q^) 
+C' i q 2 q' 2 + C' & q\' 2 + C'iq 2 q" + ---\ 

(6) 

The counter term coefficients will be determined by stan- 
dard x-square fitting procedure so that the difference be- 
tween Viow-k and (Vbare+Vcounter) is minimized. 

We note that for I and I' being both S-wave, the leading 
term for V coun t e r is Co which is momentum independent 
and corresponds to a delta interaction. For other chan- 
nels, however, the leading term for the counter potential 
is momentum dependent. For example, for I and V being 
both P-wave, the leading term is qq'Co. 

If we take the projection operator P to project onto 
states with momentum less than A, namely P = X)/c<A I 
k)(k |, we can rewrite Vbare of Eq.(5) as PVbareP- From 
now on we shall also abbreviate V coun t e r as Vct- The 
above counter term approach for Vi ow -k is similar to 
that used in RG-EFT 0, 01, as both methods impose 
some cutoff, integrate out the modes above that cutoff, 
and then attempt to recover the information contained 
in those states in a counter term series. Ultimately, the 
main difference is the the cutoff scale used and the result- 
ing potential to which Vct is added. As already noted, 
in the representation of Vi ow -k given above, we use a mo- 
mentum cutoff above which no constraining data exists, 
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and thus our "slow" potential is a bare NN potential 
projected onto the low momentum states. In the EFT 
scheme, the "slow" potential is obtained by integrating 
out heavy dynamical degrees of freedom above the chi- 
ral symmetry breaking scale. Hence, counter terms are 
added to Vu g ht, which is given by light mesons below 
the cutoff scale, and usually taken to be V one -pion with 
unrestricted pion momentum. 

We have also tried to fit Vi ow -k solely to a low order 
momentum expansion, and found the results to be less 
than satisfactory; the number of terms needed to produce 
an accurate representation is much higher than the order 
used here. Attempts at fitting (V one -pion+ Vct) also 
generated results far less satisfactory than fitting it by 
(Vbare+Vcr)- We shall discuss these points later. 

As mentioned earlier, the counter term coefficients are 
determined using standard fitting techniques. We per- 
form this fitting over all partial wave channels, and find 
consistently very good agreement. In Fig. 1 we compare 
some 1 S'o and 3 Si matrix elements of (PVbareP + Vct) 
with those of Vi ow -k for momenta below the cutoff A. A 
similar comparison for the 3 Si — 3 Di and 3 Z?i channels is 
displayed in Fig. 2. We have also obtained very good 
agreement for the P-waves, as displayed in Fig. 3. It 
may be mentioned that here the momentum factor qq' 
of Eq.(6) is essential in achieving the good fit as shown. 
As demonstrated by these figures, the two methods yield 
nearly identical effective interactions, lending strong sup- 
port that Viow-k can be very accurately represented by 
(PVbareP +Vct)- Furthermore Vct is a very short ranged 
effective interaction. 

To further check the accuracy of our counter term 
approach, we have calculated the phase shifts given by 



Viow-k and compared them with those given by (Vbare + 
Vct), as seen in Fig 4, where we have found that the 
phase shifts are almost exactly preserved by the counter 
term approach. The phase shifts are plotted up to a lab 
energy of 325 MeV, and only at this high momentum, 
do they begin to differ slightly from those obtained from 
V ow -k- In passing, we mention that the deuteron binding 
energy is also accurately reproduced by our counter term 
approach (e.g., for the CD Bonn potential the deuteron 
binding energy given by the full potential versus the 
counter term approach are respectively 2.225 and 2.224 
MeV). Thus, we can accurately reproduce Vi ow -k and all 
its results by using {PVbareP + Vct)- 

Now let us examine the counter terms themselves. In 
Table I, we list some of the counter term coefficients, us- 
ing CD-Bonn as our bare potential. It is seen that the 
counter terms are significant only for the S, P and D par- 
tial waves. We do not list results beyond the Z P^ partial 
waves, as the coefficients for them are all zero up to the 
level of our numerical accuracy (all entries in the table 
with magnitudes less than 10 -4 have been set to zero). 
It is of interest that except for the above partial waves, 
Viow-k is essentially the same as PVbareP alone. This be- 
havior is clearly a reflection that Vct is basically a very 
short ranged effective interaction. From the table, Co is 
clearly the dominant term in the expansion. Coefficients 
beyond C4 are generally small and can be ignored. In the 
last row of the table, we list the rms deviations between 
Vi ow -k and PVbareP + Vct', the fit is indeed very good. 

Comparing counter term coefficients for different Vbare 
potentials can illustrate key differences between those po- 
tentials. Thus, in Table II, we compare the low order 
counter terms obtained for the CD-Bonn |2l|, Nijmegen 



5 



TABLE I: Listing of counter terms for all partial waves obtained from the CD-Bonn potential using A = 2fm _1 . The units for 



the combined quantity q l q n C„q n of Eq.(6) is fm, with momentum q in fm 1 . All counter terms for higher order partial waves 
are zero up to our numerical accuracy. 



Wave 



Co 



C 2 



G'a 



c:; 



A r , 



'So 
3 S 1 



-0.1580 
-0.4651 



4 Si - d Di 2.879E-2 



-1.943E-3 



i D 1 
'Pi 



Pi 

>P2 



-1.309E-2 
5.884E-2 
-3.581E-3 
-1.918E-4 



-4.311E-2 -7.305E-4 
-5.566E-2 4.821E-4 
-5.479E-2 -5.916E-4 
-1.294E-2 3.075E-4 



3.561E-4 
-1.824E-3 
1.969E-3 
1.709E-4 
7.999E-4 
2.874E-4 
8.636E-4 
7.469E-4 



■8.469E-4 -1.191E-4 
-1.163E-2 -3.243E-4 5.181E-4 
-3.676E-3 -1.376E-4 6.300E-4 











-2.049E-4 -2.594E-4 



1.098E-4 







-2.065E-4 -2.721E-4 
-4.289E-4 -2.179E-4 



-1.191E-4 
5.181E-4 
-2.449E-4 













0.0002 
0.0003 
0.0025 

<0.0001 
0.0002 

<0.0001 
0.0002 
0.0002 



TABLE II: Comparison of counter terms for Vi ow -k obtained 
from the CD-Bonn, Argonne V18, Nijmegen, and Paris poten- 
tials using A = 2fm _1 . The units for the counter coefficients 
are the same as in Table I. 



Co 


c 2 


c 4 






'So -0.158 


-0.0131 


0.0004 


-0.0008 


CDB 


-0.570 


0.0111 


-0.0005 


0.0004 


V18 


-0.753 


-0.0099 


0.0003 


0.0002 


NIJ 


-1.162 


-0.0187 


0.0004 


-0.0002 


PAR 


3 Si -0.465 


0.0588 


-0.0018 


-0.0116 




-1.081 


0.0822 


-0.0002 


-0.0107 




-1.147 


0.0682 


0.0004 


-0.0100 




-2.228 


0.0251 


0.0013 


-0.0103 




3 Si- 3 £>i 0.0288 


-0.0036 


0.0020 


-0.0037 




0.0209 


-0.0025 


0.0018 


-0.0037 




0.0242 


-0.0006 


0.0040 


-0.0094 




0.0184 


-0.0015 


0.0019 


-0.0061 





|2^. Argonne H3 and Paris 0] NN potentials. The C 
coefficients for these potentials are significantly different, 
indicating that one of the chief differences between these 
potentials is the way in which they treat the short range 
repulsion. For instance, the Paris potential effectively 
has a very strong short-range repulsion and consequently 
its Co is of much larger magnitude than the others. 

Our counter term Vct has been determined by re- 
quiring a best fit between (Vb are +VcT) and Vfow-fc- 
We have explored other schemes of fitting: We have 
tried to determine Vct by requiring a best fit between 
(Vone-pion+Vcr) and Viow-k- But the results are far 
from satisfactory the resulting rms deviation at best fit 
being too large. It may be of interest to determine Vct 
by requiring a best fit between Vi ow -k and Vct alone. 



We have also tried this, with similar unsatisfactory re- 
sults (large resulting rms deviation). A possible expla- 
nation may be the following. A main portion of Vi ow —k 
come from Vt wo - P ion and other higher order processes; it 
appears that such contributions can not be compensated 
by the Vct of the simple low-order form of Eq.(6). 



IV. SCALING OF COUNTER TERMS 

In Fig. 5 we display the scaling behavior of counter 
terms Co and C2 with respect to the decimation mo- 
mentum A. We note that for the 1 5'o channel, Co and 
C2 display a weak A dependence. This is a welcome re- 
sult. We should note that in MEEFT, we are obliged to 
choose A = 2.1fm _1 , because phase shifts from experi- 
ments have been carried out up through this momentum 
and the more effective in MEEFT means that we fit all 
available experimental data by Vi ow -k- However, the Co 
for the 3 Si channel varies significantly with A. In this 
channel, there is tensor force which is a mid-range inter- 
action coming from 7r and p mesons; it has large momen- 
tum components in the intermediate momentum region 
of several fm -1 . As we lower A through this region, we 
are actually integrating out a predominant portion of the 
tensor force, in addition to the short range repulsion of 
Vbare- As a result, we would expect the Co to change 
to compensate for this loss of the tensor force. In con- 
trast, the Co for 1 Sq comes mainly from integrating out 
only the short range repulsion, and little variation is seen. 
Conversely, as A increases, the tensor force is largely re- 
tained within the cutoff, and the two Co's come close 
to each other as seen in the figure. To justify this line 
of thinking, we removed the tensor force to see how it 
would effect the 3 Si Cq. We found that the scaling for 
this term was indeed quite different, displaying a flat be- 
havior over the entire A range; Co only changed from 0.05 
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FIG. 5: Scaling behavior of S-wave leading order counter 
terms. 

at A = 1.2fm to 0.07 at A = 3.0fm. 

It should be observed that our Fig. 5 does not cover 
the small A region. The reason for this is that we have 
found that our counter term approach is not applicable 
to cases with A ~ 1/m -1 or smaller. In this region, we 
have not been able to achieve a satisfactory fit between 
Viow~k and (Vbare + Vcr)- A likely reason for this is that 
for small A's, we are integrating out a major part of the 
tensor force, which can not be compensated by the simple 
low order counter terms as indicated by Eq.(6). 

As we mentioned earlier, the counter terms are all 
rather small except for the S waves. This is consistent 
with the RG-EFT approach where the counter term po- 
tential is mainly a delta function force When using 
A ~ 2/m~ 1 , we see that the Cq coefficients for the 1 Sq 
and 3 S\ channels are different. This suggests that the 
counter term potential is spin dependent and may be 
written as 

V CT « (/o + f a ■ a)S(r), (7) 

with f Q = (3C ( 3 Si) + C (%))/4 and ft = (C ( 3 ^) - 
Co( 1 S'o))/4. For the coefficients of Table I, we have ft = 
-0.382 fm and ft = -0.077 fm. As indicated in Fig. 5, 
this spin-dependent factor ft will diminish as A increases. 



V. COMPARISON OF SHELL MODEL MATRIX 
ELEMENTS 

In shell model calculations, a basic input is the rela- 
tive matrix elements (nlSJT \ V e ff \ n'l'SJT), where 
nl and n'V denote the oscillator wave functions for the 
relative motion of the two nucleons and S and T are the 



TABLE III: Shell model relative matrix elements of Vi ow -k 
and (PVnnP + Vct) (column CT), with A = 2.0 fm" 1 and 
TiLu =19 MeV. Matrix elements are in units of MeV. 



TSj nl 


n'V 


CDB 


CDB-CT Nij 


V18 


100 


00 


00 


-9.053 


-9.053 


-9.066 


-9.062 




10 


00 


-6.744 


-6.742 


-6.698 


-6.609 




10 


20 


-3.527 


-3.529 


-3.363 


-3.375 


ni i 

Ull 


nn 
uu 


00 


-11.00 


1 1 nn 

-11. uu 


i n 79 

-1U. i A 


-10.80 




02 


no 

uu 


-8.440 


-8.405 


-8.636 






02 


02 


2.450 


2.453 


2.253 


2.263 




10 


20 


-6.873 


-6.850 


-7.033 


-6.738 




10 


22 


-6.969 


-6.940 


-7.562 


-7.637 




12 


22 


2.281 


2.278 


2.070 


2.079 


001 


01 


01 


3.835 


3.858 


3.587 


3.670 




11 


01 


3.757 


3.781 


3.446 


3.440 


110 


01 


01 


-2.482 


-2.452 


-2.407 


-2.557 




11 


01 


-1.124 


-1.099 


-0.900 


-1.162 



two-nucleon spin and isospin. The relative momenta 1(1') 
and S couple to total angular momentum J. From these 
matrix elements, the two-particle shell model matrix el- 
ements in the laboratory frame can be calculated [38| . 
Thus if Vi ow -k and PVb are P + Vct are to prove useful in 
a broad sense, they must preserve these matrix elements. 

In Table III we provide a comparison between these 
matrix elements as obtained from several input poten- 
tials. To begin, we note that the harmonic oscillator 
matrix elements for Vi ow -k are approximately the same 
regardless of which of the four bare potentials we use. 
So, as far as providing input for shell model calculations, 
this is further evidence that the Vi ow -k is approximately 
unique. Next, we come to the main result of the ta- 
ble, which is that the matrix elements for Vi ow -k and 
Vbare + Vct are in very good agreement. We have used 
the CD Bonn for the bare potential, but this holds for 
the other bare potentials as well. In fact for the 1 So, 
3 Si, and 3 Di channels, the elements are virtually exact 
out to four significant figures, and there is only minimal 
disagreement in the coupled channel. 



VI. CONCLUSION 

While there are several models for the nucleon-nucleon 
potential in use today, they suffer from uncertainty and 
model dependence. Motivated by RG-EFT ideas, a low 
momentum NN interaction, Vi ow -k, has been constructed 
via integrating out the high momentum, model depen- 
dent components of these different potentials. The result 
appears to give an approximately unique representation 
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of the low momentum NN potential, and the main issue 
addressed in this paper was whether Vi ow -k could accu- 
rately be cast in a form V\, are + Vct-, where Vct is a low 
order counter term series, and what the physical signif- 
icance of this counter term series was. We have shown 
that this was indeed the case as Vi ow -k is nearly identical 
to Vnn+Vct over all partial waves, and that Vnn + Vct 
reproduces both the deuteron binding energy and the 
NN phase shifts in the low momentum region. We have 
found that only the leading terms in the counter term se- 
ries have much significance, indicating that the counter 
term potential is mainly a short range effective inter- 
action which can be accurately represented by a simple 
low-order momentum expansion. Furthermore, we exam- 
ined the scaling properties of the S-wave counter terms 
with respect to A and found that the tensor force is re- 
sponsible for the behavioral differences exhibited between 
the 3 Si channel and the 1 Sq channel. Finally, we exam- 



ined the potential for using Vi ow -k and Vnn + Vct in 
shell model calculations. Again, it was seen that not only 
does Vio W -k give approximately identical input matrix el- 
ements irrespective of which bare potential was used, but 
that the counter term approach provided nearly the same 
results. Thus we have shown that the high momentum 
information integrated out of each bare NN potential can 
be accurately replaced by the counter terms, making the 
use of V ow -k in a broad context much more tractable. 
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